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ABSTRACT 

Context. EX Hya is one of the few double-lined eclipsing cataclysmic variables that allow an accurate measurement of the binary 
masses. 

Aims. We analyze orbital phase-resolved UVES/VLT high resolution (A/AA 27000) spectroscopic observations of EX Hya with 
the aims of deriving the binary masses and obtaining a tomographic image of the illuminated secondary star. 

Methods. We present a novel method for determining the binary parameters by directly fitting an emission model of the illuminated 
secondary star to the phase-resolved line profiles of NaLt8183/8195 in absorption and emission and CaII/t8498 in emission. 
Results. The fit to the Nal and Call line profiles, combined with the published K[ , yields a white-dwarf mass M\ = 0.790 ± 0.026 M Q , 
a secondary mass M 2 = 0. 108 ± 0.008 M Q , and a velocity amplitude of the secondary star K 2 = 432.4 ±4.8 km s~ ] . The secondary is 
of spectral type dM5.5 ± 0.5 and has an absolute A"-band magnitude of Mk = 8.8. Its Roche radius places it on or very close to the 
main sequence of low-mass stars. It differs from a main sequence star by its illuminated hemisphere that faces the white dwarf. The 
secondary star contributes only 5% to the observed spin-phase averaged flux at 7500A, 7.5% at 8200A, and 37% in the A^-band. We 
present images of the secondary star in the light of the Nal doublet and the Call emission line derived with a simplified version of 
Roche tomography. Line emission is restricted to the illuminated part of the star, but its distribution differs from that of the incident 
energy flux. 

Conclusions. We have discovered narrow spectral lines from the secondary star in EX Hya that delineate its orbital motion and allow 
us to derive accurate masses of both components. The primary mass significantly exceeds recently published values. The secondary is 
a low-mass main sequence star that displays a rich emission line spectrum on its illuminated side, but lacks chromospheric emission 
on its dark side. 

Key words. Methods: data analysis - Stars: cataclysmic variables - Stars: atmospheres - Stars: fundamental parameters (masses) - 
Stars: individual (EX Hya) - Stars: late-type - Stars: white dwarfs 



1. Introduction 

EX Hydrae (orbital period P or b = 98min) is the prototype of 
the short period version of intermediate polars, the subclass 
of cataclysmic variables in which a magnetic white dwarf ac- 
cretes from a surrounding gaseous disk or ring. Absorption 
or emission lines from the secondary star have not been con- 
vincingly detected so far, because of the strong veiling from 



the accretion disk a n d magnetic funne l dPh lion et al.l Il997t 



lEisenbartetail 120021: IVande Putte etal.l 120031) . Consequently, 
repo rts on the masses of the components have been controver- 
sial dHellier et alJ|1987t iFuiimoto & Ishidalll997t iHurwitz et al. l 



1997 



2003 



Beue rmann etalj |2003: Hoo gerwerfet al. 2004). 



In this paper, we report high-resolution phase-resolved 
blue and red spectrophotometry of EX Hya that reveals nar- 
row lines from the secondary star, notably KI/17665/7699 and 
NaI/18183/8195 in absorption and emission, a narrow emission 
component of Call/18498, and a forest of faint emission lines 
from neutral and singly ionized metals. These lines combine to 
define a unique value for the radial velocity amplitude of the 
secondary and, combined with the published K\ , allow us to de- 



rive accurate masses of both binary components. Our approach 
involves a Roche tomographic analysis of the illuminated sec- 
ondary star and the synthesis of the complex line profiles. The 
primary mass of 0.79 M Q implies a mass transfer rate close to 
that expected from gravitational radiation. 

The broad emission lines from the accretion disk and funnel 
that dominate the blue spectra will be discussed elsewhere. 

2. Observations 

EX Hya was observed in service mode with the UVES spec- 
trograph at the Kueyen (UT2) unit of the ESO Very Large 
Telescope, Paranal/Chile, on January 23 and 26, 2004. Table 1 



Table 1. Journal of the VLT/UVES observations of EX Hya. Red 
and blue spectra were acquired simultaneously (see text). 



* Based on observations collected with the ESO Very Large 
Telescope, Paranal, Chile, in program 072.D-0621(A). 



Date 


Target 


UT 


Number of 


Exposure 








spectra 


(s) 


Jan 23, 2004 


EX Hya 


5:42-8:17 


26 


300 


Jan 26, 2004 


EX Hya 


6:25-8:38 


22 


300 


Jan 29, 2004 


G1300 


11:09-11:25 


4 


60/300 
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Fig. 1. Flux cal ibrated spin-averaged b lue and red UVES spectra of EX Hya (black curves) superp osed on the overa ll spectral energy 
distribution of Eisenbar t et al.1 (l2002|) (green curves) . Also shown is the spectrophotometry of iBath et al. I (1 198 lh (small dots, see 



text), the Walraven photometry of (ISiegel et al.l l 1989) (open circles), and unp ublished UBVRIJHK photometry (open triangles, see 
text). Also shown is the optical/infrared spectrum of the dM5.5 star G1473AB (Eisenbart et al. 2002) adjusted to the flux level of the 
secondary star in EX Hya. 



contains a log of the observations. We adopted the ESO UVES 
pipeline reduction that provides fiat fielding, performs a num- 
ber of standard corrections, extracts the individual echelle or- 
ders, and collects them into a combined spectrum. The flux cal- 
ibration, except for the seeing correction (Sect. 13.21 1. is also part 
of the pipeline reduction. The wavelength calibration is derived 
from spectra taken with a ThAr lamp. Blue and red spectra were 
measured simultaneously in the wavelength ranges 3760-4980A 
and 6706-8521 A with pixel sizes of 0.030 and 0.041 A, respec- 
tively. With a slit width of 1", the FWHM resolution is 0.175A 
and the spectral resolution is A/AA ^ 47 000. To improve the 
S/N ratio, we rebinned the data into 0.3A bins (~ 10km s~') 
and obtained an effective resolution A/AA ^ 27 000. Exposure 
times were 300 s with dead times of typically 60 s between ex- 
posures due to readout. The resulting orbital phase resolution is 
A098 = 0.061, the spin phase resolution is A<p(,j = 0.090. For the 
subsequent analysis, we combined the data of both nights into 
16 orbital phase bins. 

Figure Q] (left panel) shows the flux calibrated and spin- 
averaged blue and red UVES spectra of EX Hya before the 
correction for telluric-line absorption. The spectra are over- 
laid on the spin-av eraged spectral energy distrib ution from 
lEisenbart et al.l d2002l) and the spectrophotometry of Bath et al~l 
( 1981). The latter is converted from spin phase 067=0.15 to 0.25 
(i.e., to spin average) using the known wavelength dependence 
of the spin modulation. Also shown is the Walraven photome- 
try of ISiegel et ail d!989h and the spin averages of unpublished 
phase-resolved simultaneous UBVRIJHK photometry taken in 
1982 by Joachim Krautter and Nikolaus Vogt (private communi- 
cation). The flux level of the mean red UVES spectrum agrees 
with the earlier data to better than 10%, whereas a wavelength- 
dependent correction between -5% and 20% was needed to ad- 
just the mean blue UVES spectrum to the Eisenbart et al. and the 
Bath et al. spectrophotometry. The high degree of internal con- 
sistency between the various measurements provides confidence 



in the flux calibration of the red UVES spectra discussed in this 
paper. 

Spectra of the M4.25 dwarf G1300 were taken with the same 
setup on January 29. A UVES spectrum of the M6 dwarf G1406, 
more akin to the secondary star in EX Hya than G1300, was 
kindly provided by Ansgar Reiners. Flux calibrated red spec- 
tra of these M-stars (corrected for telluric lines) are displayed in 
Fig. |2] The gap in the G1406 UVES spectrum between 8190 and 
8403A is filled in with a medium-resolution archival spectrum 
of ours (dotted curve). The red spectrum of G1300 fits the Kron- 
Cousins 7 c -band photometry (I c = 9.17, Leggett 1992) without 
any adjustment confirming the accuracy of the UVES flux cal- 
ibration, whereas the blue spectrum shows a similar moderate 
loss of blue light as noted above. The G1406 spectrum is ad- 
justed to I c = 9.39 (Leggett 1992). As discussed below, the sec- 
ondary star in EX Hya is of spectral type dM5 to dM6 and con- 
tributes about 5% of the mean flux of EX Hya at 7500A. The 
correspondingly adjus ted spectrum of the dM5.5 star G1473AB 
( lEisenbart et al.ll2002h illustrates the dominance of the disk and 
funnel emission in EX Hya (Fig.Q]). 



3. Data Analysis 

In this section, we discuss the accuracy of the wavelength cali- 
bration, outline the seeing correction that is part of the flux cal- 
ibration, and discuss the correction of the individual red spectra 
for telluric line absorption. 



3.1. Wavelength calibration 

The UVES pipeline reduction provides wavelength-calibrated 
spectra. The accuracy of this wavelength calibration can be ver- 
ified in the red spectral region by measuring the positions of the 
numerous telluric lines of O2 and H2O. As an example, Fig. [3] 
shows the water vapor lines near 8200A. Comparison of theo- 
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Fig. 2. Flux calibrated red UVES spectra of the M6 dwarf G1406 
and the M425 dwarf G1300 (shifted upwards by two units). The 
'residual' spectra shown as the two bottom curves refer to our 
method of analyzing the EX Hya spectra and are explained in 
Sect. 13.31 The bottom one refers to G1406 and is shifted down- 
wards by one unit. 



retical and observed positions of unblended lines with moder- 
ate optical depths yield a mean difference in these positions of 
AA = 3 ± 3 mA throughout the red spectral range and a stan- 
dard deviation of <j — 13 mA. Hence, the wavelength calibration 
is quite accurate and the systematic error in the derived radial 
velocities does not exceed 1.0 km s -1 . 



3.2. Seeing correction 

Seeing information is provided by the DIMM on Paranea l 
(Difference Image Motion Monitor, Sarazin & Roddier 1990) 
and by the width of the individual spectral images perpendicular 
to the dispersion direction. Both methods yielded compatible es- 
timates of the FWHM seeing that varied between 0.9" and 2.2" 
on January 23 and between 0.7" and 1 .2" on January 26. We cor- 
rected each spectrum for the seeing losses assuming a Gaussian 
point spread function and a central position of the source in the 
1 .0 arcsec slit. This correction is part of our flux calibration and 
is included whenever we quote absolute fluxes. In particular, this 
holds for the spectra in Fig. [TJ and for the light curves shown in 
Fig- Eland discussed in Sect. 14.41 below. 

3.3. Telluric line correction 

We reconstruct the stellar flux incident on the atmosphere by 
correcting the observed spectra for the absorption in the telluric 
lines. Given the structure of the spectrum as shown in Fig. [3] this 
is a formidable task, in particular, since the water vapor content 
of the atmosphere on January 23 and 26, 2004, was rather high. 
Absorption in the strongest H2O lines near 8200A reached 88%, 
whereas the depth of the NaL18183/8195 doublet is only 2.6% 
of the continuum. 

Telluric correction of the observed spectrum Fi{A) of orbital 
phase bin i requires division by exp(-Ti(/l)), where T;(/t) is the 
wavelength dependent optical depth of the atmosphere. We need 
to define a template t(A) that allows us to construct the individ- 
ual r;(/l) and to reach the desired accuracy of the reconstructed 
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Fig. 3. Mean normalized spectrum of EX Hya in the vicinity of 
the NaI/18 183/8 195 doublet. Wavelength bins are 41 mA and the 
FWHM spectral resolution is 175 mA. The two spectra shown 
at the top are explained in Sect. 13.31 The Nal doublet becomes 
detectable only after telluric correction. 



fluxes. To start with, we employed the normalized spectrum of a 
featureless standard star obtained in a different night as template 
t(A). This proved unsuccessful, because the optical depth in wa- 
ter vapor was only about 1/3 of that in our observations and the 
optical depth structures did not match. We resorted, therefore, 
to an entirely different method using the orbital mean red spec- 
trum of each night as a telluric template. At first glance, this 
might seem inappropriate, because it removes all spectral struc- 
ture from the mean of our set of phase-resolved spectra. For the 
special case of a cataclysmic variable with its rapid variability, 
the method is attractive, however, because it removes all spectral 
structure that does not display orbital (or more generally tempo- 
ral) variability. Although the mean spectrum vanishes, the indi- 
vidual corrected exposures that we refer to as 'residual' spectra 
retain much of the spectral structure at wavelength scales less 
than that of the radial velocity variation of ~ 20A in case of EX 
Hya. 

To explain the method, we present a model calculation with 
the 4- A binned spectra of G1300 and G1406 shown in Fig. [2] Let 
F(A) be the spectrum of the M-star and F[(A) the same spectrum 
for orbital phase bin i shifted in wavelength according to the ra- 
dial velocity of the secondary star in EX Hya. The mean of N 
phase bins is F(A) = (l/N) F[(A) and the residual spectra are 
fi(A) = F\(A) - F(A). Figure|2](bottom curves) shows the residu- 
als of G1300 and G1406 at orbital phase 98 ,i = 0.75, calculated 
for an orbital velocity amplitude of 432.4km s _! and a spectral 
flux that does not vary with phase. Spectral structure that ex- 
tends over AA > 20A is lost, while structure on a shorter scale is 
preserved with an amplitude up to 90% of the original. 

For EX Hya, we express the 'raw' spectrum Fi(A) of orbital 
phase bin 2 as 



F i ^) = C i U)n i (A), 



(1) 



where Ci( A) represents the smooth wavelength dependence of 
the continuum and ti{{A) is the normalized spectrum that contains 
all spectral structure. We use the orbital mean of the n\(A) as the 
telluric template n(A) = exp(-r(/l)). The correction function for 
phase bin i is calculated as exp(-T;(/l)) = exp(-ajr(/l)) = n(A) a ' 



4 
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Fig. 4. Flux-calibrated light curves of EX Hya obtained on 
January 23, 2004. The bands U' and B refer to the mean flux in 
the intervals 3770^000A and 4000^800A, respectively. The 
integrated flux of the HB emission line is given in units of 
10~ 13 ergcm"" 2 s -1 . The triangles indicate the t imes of spin max- 
ima p redicted by the quadratic ephemeris of iHellier & Sproatsl 
(1992). 



with a parameter a\ that is independent of wavelength. We de- 
termine the Ci(A) by fitting a low-order polynomial to the con- 
tinuum of the raw spectra. The template is then fitted to the re- 
sulting rii(A) over a specified spectral range allowing for a small 
wavelength shift 5A-, (of order mA). Hence, the residual spectrum 
for phase bin i is 

MA) = Fi(A)[n(A + 6A i )r a> - Q(A) , (2) 

with fit parameters a\ and 5A[. In our data, the optical depth 
varies little with orbital phase and the a\ stay close to unity. 
We define the continua C;(/l) by quadratics fitted to the F\(A) 
at 6800A, 7720A, and 8362 A for A < 8000 A and at 7130A, 
7854A, and 8362A for A > 8000A. At these wavelengths, any 
temporal variation of the spectral flux is removed and the resid- 
ual spectra equal zero. Any remaining orbital modulation of the 
residual flux is relative to these wavelengths. Telluric absorption 
is considered separately for absorption by O2 and H2O lines. 
The template is fitted to the O2 lines between 7620A and 7657A 
avoiding the KI lines and the deepest part of the A-band at 7593- 
7617A. It is fitted to the H 2 lines over the bands 8000-8150A 
and 8220-8350A avoiding the Nal doublet. We restrict the sub- 
sequent analysis to A > 7300A and piece the residual spectra 
together from the 02-corrected part below 8000A and the H2O- 
corrected part above 8000A. The two fits remove the weaker tel- 
luric lines outside the immediate fit intervals almost perfectly 
and only the 7593-7617A section remains problematic, where 
absorption at the line centers reaches 99.8%. The spectra of both 
nights are then combined and collected into 16 phase bins. The 
result is shown in Fig. [5] Spectral structure on a short wave- 
length scale is preserved in the individual reconstructed spectra 
although they average to zero. 



3.4. Statistical noise 

A natural consequence of any telluric-line correction is the in- 
creased noise level of the restored fluxes at the positions of the 
strong lines. We opted to keep all data points in the resulting 
spectra and account for the enhanced noise at the telluric line 
positions by the increased statistical errors of the affected data 
points. The orbital mean of the reconstructed normalized spec- 
tra for the first night is shown in Fig. [3] (second curve from top, 
shifted upwards by 0.2 units). As expected, it equals unity with 
only minute deviations that disappear, too, when it is rebinned 
to 0.3A. The statistical error of the individual 0.3A spectral bins 
is determined from the mean rms noise in regions that avoid the 
strongest telluric lines and amounts to 0.7% of the continuum 
(Fig. [3] uppermost curve). It reaches 2% in the deepest H2O lines 
around 8200A and huge values in the deepest lines of the A-band 
between 6593 and 7617A. The 0.7% level corresponds to an ab- 
solute error of 6 x 10~ 17 erg cm _2 s _1 A -1 near 8200A. 

4. Observational results 

The gray plots of Fig. [5] show the 16 phased-resolved resid- 
ual spectra for the wavelength intervals 4875^975A and 
7300-8520A displayed twice. The wavelength bins are 0.3A. 
In the selected spectral ranges, the strongest broad emis- 
sion lines from the accretion disk and funnel are HeI/14922, 
OLT7772/7774/7775, Paschen P17, CaH/18498, and Paschen 
P16. The narrow spectral features from the secondary star 
include the TiO band heads at 7589A and 8432A, the 
KLT7665/7699 and the NaLl8183/8195 lines in absorption and 
emission, Call/18498 in emission, and numerous metal emis- 
sion lines, among the stronger ones FeI/14957, FeI/18327 and 
FeI/18388. The absorption lines are strong between <p^ = 0.65 
and 1.35, whereas the emission components reach peak flux at 
098 — 0.50, both indications of an origin from the illuminated 
secondary star. The narrow and strong CaII/}8498 emission line 
is superimposed on a complex background of the much broader 
Paschen P16, P17, and Call/18498 disk lines. We have separated 
the narrow line interactively from the background and included 
an estimate of the uncertainty of this procedure in the errors of 
the individual data points. 

4.1. The secondary star in EX Hya 

Figure [6] displays mean spectra for the phases when the illumi- 
nated face of the secondary star or its dark side are in view (black 
curves). The spectra are summed over the respective phase in- 
tervals, are placed on an absolute flux scale, and shifted into 
the rest systems of the emission or the absorption lines, respec- 
tively. To appreciate the weakness of the lines note that the depth 
of the Nal doublet is 2.6% of the continuum at 8200A (Fig.[TJ. 
Information on the spectral type of the secondary star in EX Hya 
can be obtained from a comparison of these spectra with those 
of the M-stars G1300 and G1406 (green curves). We have used 
the M-star spectra at their full spectral resolution, adjusted their 
fluxes with the ratio of the solid angles Q, 

F S ec = ^M^sec/^M = Fu(fisecl 'd) 2 (dM/ 'Rm) 2 '\ (3) 

formed the residuals as explained in Sect. l3.3l and rebinned them 
to 0.3A. Here, dyi and Rm are the distan ce and the stellar ra- 
dius of the M-st ar (Beuerma nn et al.lll999l) . whereas d - (64.5 ± 
1.2)pc dBeuermann et al.l2003l) and R< P r = (0.15 16 + 0.0034) R Q 
(see Sect. 16.31 below) denote the distance of EX Hya and the 
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dwarf G1406 and the M4.25 dwarf G1300. 



mean radius of its secondary, respectively. The ratio of the solid 
angles amounts to Q. se J£l M = 0.00421 + 0.00036 for G1300 and 
fi sec /Q M = 0.00170 + 0.00015 for G1406. The residual spectrum 
of the dark side of EX Hya is displayed twice in the three pan- 
els of Fig. |6j at the center and the bottom. For A < 7605 A, we 
have chosen the phase interval 098=0.76-0.91, when the emis- 
sion lines have already disappeared and the blueshift has moved 
the TiO band head at 7589A away from the telluric A-band. For 
A > 76 15 A, we sum over the entire dark side of the secondary 
(098=0.76-1.18). Our comparison indicates that the residual 
spectrum of G1406 reproduces the TiO band head of EX Hya 
at 7589A perfectly. The KI lines of G1406 fit almost as well, 
whereas those of G1300 utterly fail. The Nal doublet suggests a 
spectral type between G1300 and G1406 and the TiO band head at 
8432A prefers G1300. Many of the low-amplitude wiggles in the 
EX Hya spectrum faithfully reproduce rotationally smoothed ab- 



sorption line structures of the two M-stars. This agreement is lost 
for A > 8445 A, where the broad disk lines veil the spectral fea- 
tures of the secondary star. The TiO band head at 8432A is also 
affected by the two flanking Til emission lines and, assigning it 
a lower weight, we conclude that the EX Hya spectrum is much 
better represented by G1406 (dM6) than by G1300 (dM4.25). The 
implied spectral type of the secondary is slightly earlier than M6 
and we settle for dM5.5 ±0.5. Interpolating between the adjusted 
fluxes of G1300 and G1406, we find that the the secondary con- 
tributes 0.92+0.28 mJy and 1 .40+0.48 mJy at 7500A and 8200A, 
respectively. This corresponds to (5.3 ± 1.6)% and (7.5 ± 2.6)% 
of the spin-averaged flux of EX Hya, respectively (see Fig.Q]). 

Although the dark side of the secondary looks like a dM5.5 
star, this does not hold for the illuminated side as shown by the 
spectrum for 098= 0.40-0.63 in Fig.|6](top). This phase interval 
is characterized by the peak flux of the emission lines and a dis- 
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appearance of the TiO/17589 absorption edge. Figure|7]shows the 
light curves of the integrated flux of the Call/18498 emission line 
(upper panel) and the difference AF 7550 of the residual fluxes at 
7550 and 7720A (lower panel). The minimum of the latter coin- 
cides with the emission line maximum stressing the weakness of 
the TiO/17589 band head on the illuminated face, a result that is 
reminiscent of Wade & Home's (1988) finding for Z Cha. 

Applying Eq. 0) to the A"-band, allows us to estimate the 
contribution of the secondary to the infrared flux of EX Hya. The 
apparent A"-band magnitudes of G1300 and G1406, K = 6.68 and 
6.08, imply A" sec = 12.62 and 13.00 for a secondary of spectral 
type dM4.25 or M6, respectively. For the preferred spectral type 
dM5.5 + 0.5, EX Hya B has K sec = 12.89 ± 0.1 1. We can alter- 
natively s tart from the surfac e brightness calibration of field M- 
dwarfs by Beuermann ( 2006) that implies S k = 4.7 1+0.11 for a 
spectral type dM5.5± 0.5. With fl sec = (0. 15 16 + 0.0034) R Q and 
d = (64.5±1.2)pc, we find/T = S K -5 \og[(R sec /R@)(l0pc/d)] = 
12.86 ± 0.13, nearly identical with the previous number. The 
absolute magnitude of the secondary in EX Hya is Mk = 
8.83 ± 0.12, in perfect agreement with that of a main sequence 
dM5.5 star. The spin-averaged apparent A"-band magnitude of 
EX Hya is K — 1 1 .8 1 (Fig.[TJ to which the secon dary contributes 
37.6 + 4.2)%, somewhat less than estimated bv lEisenbart et al.l 
2002). These authors noted already that part of the modulation 
of the infrared bands at one half the orbital period that looks 
like the ellipsoidal modulation of the secondary star may be 
produced by the bulge on the accretion disk. We confirm this 
suspicion by noting that the optical bands of the unpublished 
UB VRIJHK photometry (Sect. |2]i display a 49-min modulation, 
too. 
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Fig. 8. Radial velocity curves of NaLi8183/8195 in absorption 
(O) and emission (+) (left hand scale) and of the CaILi8498 
emission fine (+) (right hand scale) vs. orbital phase i/?y8HS 
(lHeliier&Sproatsl fl992i. Solid lines represent sinusoids least 
squares-fitted to the data (see Sect. 14.2b . The error bars are 
smaller than the data points. 



4.2. Radial velocity amplitudes and systemic velocity 

The brightest emission lines suited for radial velocity measure- 
ments are CaILi8498, KU7665/7699, FeU8327, and FeL14958. 
Of the absorption lines, only the NaLl8183/8195 doublet yields 
reliable results. In KLJ7665/7699, the emission component is 
relatively stronger and the profiles are affected by residuals of 
the telluric O2 lines. In this initial analysis, we measured ra- 
dial velocities using Gaussians for the single lines and a dou- 
ble Gaussian with a separation of 11.57A for the Nal doublet. 
Figure [8] shows the radial velocity curves of the Call emis- 
sion line measured from the 48 individual exposures and of the 
Nal doublet measured from the 16 phase-binned spectra. The 
abscissa 098. HS is the orbital phase c alculated from the linear 
ephemeris of Hellier & Sproatsl (ll992l) . The Nal absorption com- 
ponent is disturbed by emission on the descending branch and 
becomes undetectable around <pgs — 0.5. We have measured the 
velocity amplitudes, the phase shifts, and the systemic velocities 
from sinusoidals fitted to all data points for the emission lines 
and to the data points between 098,hs = 0.71 and 1.21 for Nal. 
The fit results are listed in Tab. [2] Although these fits do not ac- 
count for the slight ellipticity visible in the data, they yield quite 
accurate values of the velocity amplitudes K' 2 of the emission 
lines and of the systemic velocity y. Neither the emission lines 
with an average K' 2 = 351kms _1 nor the Nal absorption lines 
with K' 2 ' = 460 km s -1 provide a reliable measure of the veloc- 
ity amplitude K2 of the secondary star, which we expect to fall 
in between. Note that the TiO/17589 band head is not suited to 
track the motion of the secondary star, because it is disturbed by 
the A-band around 098=0.25 and by FeI/17584/7586 emission 
between 098=0.2 and 0.8. 

The systemic velocities of the individual lines in Tab. [2] are 
generally in good agreement, although the Nal absorption line 
value may be affected by the emission component and the first 
KI value by possible blending with Fel at 7664. 29A. Excluding 
these, the mean systemic velocity of the remaining four lines 
becomes y- -59.0 + 0.9 kms _I . The small systematic error 
of our wavelength calibration adds only 0.5kms _1 to the error 
budget of y (see Sect. 13. It . Confirmation on the common 
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Table 2. Results for sinusoidal fits u ra( j = K sin[27r(</>98 + <po)] + 
y to the radial velocities of narrow emission lines and the Nal 
absorption line doublet. The errors in brackets refer to the last 
digit. See text for systematic errors. 
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-61 ± 


4 


cm 


Fel 


8327.05 
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5 


0.019(3) 
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5 


cm 


Call 


8498.02 


357 + 


3 


0.020(4) 


-57 ± 


5 


cm 


Nal 


8183.26/8194.82 


460 ± 


5 


0.021(2) 


-68 ± 


3 


abs 



origin of all narrow emission lines and on the value of y is 
obtained from the numerous lines seen in the blue spectra. 
They originate from neutral or singly ionized species and are 
restricted to the phase interval when the illuminated face of the 
secondary is in view. The better defined ones with equivalent 
widths between 0.04A and 0.11 A are HL13770, 3797,4102, 
FeLt3856, 3889, 4198, 4271 , 43 14, 4426, 4890, 4957, Call K, 
MnU4034, ArU4045, OIU4144, NaL14216, NeII/AlII/14227, 
OL14232, Neli4274, CrU4583, and MgIU4938. Narrow 
emission components are not detected in the Hel lines and 
in Hell/14886. We have studied these faint lines in a cursory 
manner shifting them into their rest system with an appropriate 
KL We obtain straight vertical lines in the 2-D image for the 
K' 2 value of the red lines, 351kms _1 . The systemic velocities 
are measured from a spectrum co-added over the phase interval 
of best visibility of the lines. Their distribution has a standard 
deviation of cr = 2.3kms~ I and a mean y = -59.6 km s~', 
indistinguishable from that of the red lines. 

Unfortunately, there are no accurate previous measures of y 
from optical lines to comp are our result with, except the rather 
uncertain values quoted by Hellie ret al.l (Il987l) . It is potentially 
important, however, to note the more positive y velocities de- 
rived from UV and X-ray emission lines that originate near the 
white dwarf surface, y = 9. 5 ± 3.0 km s" 1 (iBelle et al.l 120051) 
and y = -2.8 ± 2.3 km s -1 dHoogerwerf et all I20041T The er- 
rors quoted by these authors are statistical ones and additional 
systematic errors may affect the numbers. Nevertheless, it is in- 
teresting to note that the difference between the optical and X- 
ray/UV results is close to the gravitational redshift expected for 
the white dwarf in EX Hya (see Sect. [7). 

4.3. Check on the origin of the narrow emission lines 

IBelle et all (120031) and iHoogerwerf et al.l (|2004) have reported 
accurate values of the radial velocity amplitude of the white 
dwarf, K\ = 59.6+2.6kms 1 and^i = 58.2±3.7 km s" 1 , respec- 
tively. We show here that the fully resolved profiles of the narrow 
optical lines in our spectra independently yield a very similar, al- 
though not quite as accurate value of K\ . To this end, note that 
the measurement of two velocity amplitudes that can be assigned 
to two points on the line connecting the two stars is equivalent to 
the measurement of the radial velocity amplitudes of both stars: 
we choose the back of the Roche lobe and the Li point. The 
peak velocity in the Nal absorption line profile at quadrature cor- 
responds to the velocity amplitude of the back of the star if we 
account for some widening by pressure broadening. We choose a 
point two pixels (20 km s ) down into the absorption line profile 



and obtain K\,. i( ± = 550+ 10km s _1 . The minimum velocity in the 
Call emission line profiles at quadrature is taken to represent the 
Li point. Since the emission lines display no additional broad- 
ening, we choose the first pixel in the Call line profiles with a 
non-vanishing flux and obtain Kn = 280 ± lOkms -1 . The ratio 
of these velocities equals /^back/^Li = -^back/^Li = 1-96 + 0.08, 
with ;tback and xli the coordinates along the x-axis of the two 
extremal points on the stellar surface as measured from the cen- 
ter of gravity in units of the binary separation. Roche geometry 
then implies q = M 2 /M l = 0.134 + 0.019,x back = 1.129 + 0.005, 

K 2 = #back/*back/(l + ?) = 430 ± 9kms-\ Ki = qK 2 = 
58 + 8 km s -1 , and K x + K 2 = 488 ± 9 km s" 1 . With the orbital pe- 
riod P or b=5895.4s and the inclination i = 77.8°, the binary sep- 
aration is a = (Ki + K 2 )P mb /(2nsini) = (4.683 + 0.088)10 10 cm. 
Finally, Kepler's law gives a total mass of M - 0.88 ± 0.05 M Q 
and component masses Mj = 0.775 + 0.045 M and M 2 = 
0.104 + 0.016 M . Since /Tback and Kn were simply read from 
the observed profiles, these numbers are approximate only. The 
important point to note is the excellent agreement between our 
value of K\ and the directly measured, more accurate radial ve- 
locity amplitude of the white dwarf quoted at the top of this 
paragraph. This agreement verifies our assumption that the low- 
velocity limit in the Call line profiles at quadrature measures the 
motion of matter that is located on the secondary star near Li . 
This assumption forms the basis of the line synthesis procedure 
described in Sect. [5] below. 

4.4. Phase shifts 

All photometrically dete rmined eclipse timings re ported over the 
last decade occur near a Hel lier & Sproatsl d 19921) orbital phase 
098,hs - 0.98. The blue-to-red zero crossing of the Nal line in 
our data takes place at = 0.982±0.001 (Tab. O indicating 

that eclipse and zero cr ossing coincide wi thin about the jitter in 
the eclipse timings (e.g. Sieg el et aljl989l) . This is confirmed by 
one pronounced dip that is superimposed on the first spin max- 
imum in the night of January 26, 2004, occurs at spectroscopic 
phase zero within the error, and probably represents the eclipse 
by the secondary star (Tab. [3]). 

The spin light curves in Fig.|4]show that blue light, red light, 
and the Hfi flux display the same spin modulation within errors. 



Table 3. Barycentric timings of the spin maxima (spin phase 
zero), the eclipse center, and the zero crossing of the absorption 
lines from the secondary star (orbital phase ze ro). The O-C val- 
ues ar e given relative to the ephemerides of iHellier & Sproatsl 
dl992l) (see text). 



Band HJD+2400000 Error O-C 

(d) (d) (d) 



(a) Spin maxima 



Blue light 


53027.7670 


0.0014 


0.0113 




53027.8109 


0.0014 


0.0087 




53030.8381 


0.0014 


0.0104 


H/?flux 


53027.7640 


0.0014 


0.0083 




53027.8109 


0.0014 


0.0087 




53030.8360 


0.0014 


0.0083 


(b) Eclipse timing 








Blue light, U/3 flux 


53030.7909 


0.0014 


-0.0004 


( c) Absorption line 


zero crossing 






Nal doublet 


53030.7893 


0.0001 


-0.0012 
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The U\ B and 7500A bands are integrals over 3760-4000A, 
4000^800A, and 7450-7550A, respectively. The H/3 flux is an 
average over +40A of the line center corrected for the underly- 
ing continuum. The solid triangles on the abscissa indicate the 
times of the spin ma x ima p redicted by the quadratic ephemeris 
of iHellieF & Sproats] d!992f) . The offset between observed and 
calculated times of maxima is A067 = -0.20 + 0.03 and, in what 
follows, we use the phase conventions 



Table 4. Parameters of the irradiation model for EX Hya and 
results of the combined least-squares fit to the NaI/18183/8195 
and the Call/18498 line profiles in 16 orbital phase bins. 



"98 



LBS +0.018 



and 



»67 



t>67,HS 



0.20. 



(4) 



Irradiation of the secondary star in EX Hya varies periodically 
at the beat period between orbital and rotational periods (210 
min) and inclusion of Eq. (01) into our illumination model assures 
correct phasing. The only other post- 19 91 timings of blue ligh t 
maxima that we are aware of are those by Eisenbart et al. (2002), 
A0 67 = -0.04 ± 0.01 at HJD 2450508, and lBelle et aTT d2005), 
A0 67 = -0.12 + 0.02 near HJD 245 1685. Although substan- 
tial fluctuations occur in the timings (Hel lier & Sproatsl ll992; 
iBelle et alJl2005l) . the numbers suggest that the offset increases 
with time and that the spin-u p of the white dwarf in E X Hya is 
slightly slower than quoted bv lHellier & Sproatsl dl992l) . A more 
regular monitoring of EX Hya with a small telescope a worth- 
while undertaking. 



5. Irradiation model 

We now embark on the construction of a detailed line synthe- 
sis model for NaLI8183/8195 and CaIU8498. Our model rep- 
resents the Roche lobe by a grid of triangular surface elements, 
each of which is characterized by photospheric absorption and 
potentially by superimposed emission from a thin layer that we 
do not distinguish geometrically from the photosphere. 



5. 1 . Line profiles 

Velocity smearing is adopted as the dominant broadening mech- 
anism, supplemented by pressure broadening for the Nal absorp- 
tion lines and minimal Gaussian broadening of the Call emission 
line. Pressure broadening is represented by a Lorentz profile de- 
rived from the Nal doublet in the M4.25 dwarf G1300. The emis- 
sion lines, on the other hand, are produced in layers above the 
photosphere and probably lack significant pressure broadening. 
The adopted Gaussian broadening with a FWHM of one pixel 
(lOkrns -1 ) merely serves to smooth irregularities arising from 
the 0.01 phase bins of the model spectra. 

In principle, modeling the complex line profiles that contain 
absorption and emission components requires appropriate ra- 
diative transfer calculations for irradiated M -dwarf atmospheres 
dBrett& Smith 1 119931: iBarman et aill2004 . Results that could 
easily be implemented are not yet available, however, and we 
consider two simple cases instead: (i) the fill-up of the absorp- 
tion line with emission by adding the two components; and (ii) 
the gradual disappearance of the absorption line in the illumi- 
nated part and its replacement by emission. The difference lies 
in the absorption line wings that are retained in case (i) and prac- 
tically disappear in case (ii). Testing both models led to better 
fits with and a clear preference for case (i) (see Fig. [14] top). 
The adopted procedure is adequate for the present data, but may 
have to be replaced by a more sophisticated approach if data of 
better statistical accuracy become available. 



No. Parameter 


Symbol 


free/ 
fixed 


Value ± Error 


(a) System parameters: 








PI Orbital period (s) 


P 

1 orb 


fixed 


5895.4 


P2 Inclination 


i 


fixed 


77.8 ± 0.4 ' 


P3 Velocity amplitude (km s" 




fixed 


59.0 ± 3.0 ' 


P4 Total mass (M Q ) 


M 


free 


898+0 031 


P5 Mass ratio M 2 jM x 


q 




0.137±0.007 3) 


P6 System velocity (km s ) 


j 


free 


-58.2 ± 1.0 


P7 Phase shift 




free 


0.018 ±0.001 


P8 Phase shift 


A067 


fixed 


-0.20 ± 0.06 


P9 Ampl. of irradiation flux 




fixed 


0.15 


(b) M;//18183/8195 absorption line parameters: 




P10 Gravity darkening 


y.\ 


fixed 


0.61 


PI 1 Limb darkening 


c, d 


fixed 


see text 


( c) Call and Nal emission line paramete, 


r s: 




P12 Disk half opening angle 


P 


free 


2 1 +1 ° 


PI 3 Terminator, cos i? 


t lim 


free 


-0.06 ± 0.06 


P14 Optical depth 


T 


fixed 


> 10 


PI 5 Limb darkening 


M l,Ca 


free 


0.42 ± 0.24 4) 


P16 FWHM (km s-') 


FW 


fixed 


10 


P17 Normalized Call 


/»(#) 


free 


see Fig.rTOl 



. . . P24 emission line fluxes 

11 From FWHM of the X-ray ecl i pse for given WD radius i?[(M[). 

2) From dHoogerwerf et all2004l ; belle et ai1l2003l) . see text. 

3) If Ki is fixed in the fit, q(M, K t ) is not an independent fit parameter. 

4) Error reflects the strong correlation between iii Ca and M. 



5.2. Parameterization of the model 

Our kinematic and illumination model has 29 parameters of 
which 24 are listed in Tab. |4] and the remaining five are nor- 
malization constants explained in Sect. 15.41 Parameters that are 
kept fixed are the orbital period P 01 -b, the inclination i, the phase 
shift A</>67, the 210-min amplitude of the irradiation flux A[ n , the 
optical depth of the emitting layer on the secondary r, the in- 
trinsic FWHM of the line emission, the Nal gravity darkening 
coefficient y^, and the constants c and d in the Nal limb darken- 
ing law (all explained either above or in the next two Sections). 
We opted to provide the white dwarf r adial velocity amplitude 
Ki a 59 + 3km&- 1 (IBelle et al.ll2003t iHoogerwerf et alj|2004l) 
as a fixed input parameter and determine its influence on the er- 
rors of the fit parameters by varying it between 56 and 62 km s . 
The mass ratio q is then a derived parameter and is determined as 
q(M, Ki). Of the free parameters, the quantities M, y, and A09g 
describe system properties, whereas the disk half thickness (3, 
the limit of illumination at c\\ m , the limb darkening/brightening 
coefficient of the Call emission M^cai tne eight numerically given 
Call emissivities (or one parameter of the analytical model), and 
the five normalization constants define the emission line model 
and are discussed in Sect. 15.41 For the numerical or analytical 
versions of the emission model, 19 or 12 parameters, respec- 
tively, are fitted. 



5.3. Malabsorption 

Limb darkening coefficients are available in the literature only 
for broad photometric bands. Here, we need the limb darkening 
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yLt=COS B 

Fig. 9. Limb darkening of the angle-dependent intensity deficit 
integrated over the Nal doublet for a logg=5 star with T^g- 
3000 K (Derek Homeier, private communication). The solid line 
is a fitted square root law with fit parameters c and d given in the 
text. 



of the integrated Nal absorption line flux. We extract this infor- 
mation from the model atmosphere of an unirradiated log g-5 
star of 3000 K kindly calculated by Derek Homeier with the 
PHOENIX code. We determine the mean Nal intensity deficit 
from the angle-dependent intensity hip) by integration over the 

doublet as / Na 0") = J line Qc(m) ~ hQ*)) dX = EWQi) I c Qi), where 
fi = cos 9 with 6 the zenith angle and I c denotes the continuum 
outside the line. Figure [9] shows this quantity normalized to the 
center of the stellar disk along with a square root fit 

/NaOu)// Na (l) = l-e(l -H)~d(l ~ (5) 

dClaretl 11998b with parameters c = -0.069 and d = 1.063. 
The EW of the Nal doublet decreases slightly as the limb is ap- 
proached and varies approximately as EW(ju) = (6.72+1 .07 /_/) A. 
The EW averaged over the stellar disk is 7.4A and both, the EW 
and the effective temperature of 3000 K, are typical of a dM5.5 
star, the best estimate of the spectral type of the secondary star 
in EX Hya (see Sect.|7]and Tab.©. 

We account for monochromatic gravity darkening at 8200A 
by a linear coefficient y,\ - 0.61 (Tab. |4| that is based on a sur- 
face brightn ess vs. effect i ve te mperature relation derived from 
the results of lBeuermannl d2006l) for field M dwarfs. 

Photospheric absorption in the Nal doublet contributed by 
the /c-th surface element at orbital phase 09g,; is then given by the 
flux increment 

dF, x = - fl , dA kMk ^77 f 1 + W " 5)) ' < 6 > 

with dA k the surface area of the element, fi k = cos^ and 6 k the 
angle between the normal to the element and the line of sight at 
orbital phase 098,;, r k the radial separation of that element from 
the center of the star, and R the mean stellar radius. The flux 
increment is collected into the appropriate 0.3 A wavelength bin 
that corresponds to the radial velocity of the element as seen 
by the observer at 098,;. For a hidden element with /i k < 0, the 
flux increment is set to zero. The flux level of the Nal doublet is 
determined by the fit parameter a. d and a further parameter a rat ; 
describes the ratio of the /18183 vs. the ^18 195 lines fluxes. 



5.4. Call and Nal emission 

Parameterization of line emission from the irradiated atmo- 
sphere is more involved. We model the Call/18498 emission and 
adopt its properties for the emission of the Nal doublet with a 
scaling factor for the different emission line intensities. We as- 
sume that the angle-integrated Call emission line flux Fez of 
surface element depends on its distance r s from the irradiating 
source and on the angle of incidence § as Fc- d x f(jj)/t% with 
j] - cos In the special case of an emitted flux proportional 
to the incident energy flux, f(rf) oc 77. Model calculations for 
f{rf) in M-dwarf atmospheres are not yet available, but the case 
of Lyman line emis sion of irradiated white dwarf atmospheres 
dKonig et al. I l2006h suggests that f(j])lr\ increases rapidly as 
grazing incidence is approached (77 = 0). We have tested analyti- 
cal formulations of f(rj) and a numerical tomographic approach. 
For the latter, we fix f(rj) at /; =0 for r\ — (terminator) and 
at fio = 1 for r\ — 0.9 (L; point). In between, we define f(r\) 
by eight free parameters fz to fg at equidistant abscissa values 
77 (m) . We interpolate between the f m and regularize the distribu- 
tion of the emissivities by adding A Em=iC/m+i -fm) 2 to the^f 2 to 
be minimized, with A a Lagrange multiplier. This choice implies 
that a smoothed version of the current set of f m is used as the 
default. 

The observed light curve of the wavelength-integrated Call 
line flux (Fig. [7] upper panel) is slightly skewed with a centroid 
at orbital phase 09g = 0.5 17. We model this asymmetry by adopt- 
ing the hot spot or bulge as a second source of irradiation that we 
locate in the orbital plane at binary coordinates x = 0,y — 0.3 in 
a system with its origin at the center of gravity, unity separation 
of both stars along the x-axis. 

Three further parameters are needed to account for (i) the 
shadow cast on the secondary by the accretion disk, (ii) a pos- 
sible horizontal energy transfer across the terminator into the 
dark side of the star, and (iii) the periodic variation of the ir- 
radiation flux received by the secondary from the rotating mag- 
netic white dwarf. We collect these dependencies into three ad- 
ditional parameters: (i) the shadow has a half opening angle /3 
independent of azimuth and a sharp edge; (ii) line emission is 
allowed to extend beyond the geometric terminator for illumi- 
nation by a point source at rj — to 77 = cy lm , where c;; m is a 
small negative (or positive) quantity; and (iii) rotational modula- 
tion is modeled by a variation of the irradiation flux of the form 
F m - oc 1 + a m cos {2n 0; n ), where 0; 1T = 067 - 098 is the irradiation 
phase, maximum irradiation occurs at <p\„ = 0, and varies 
with the beat period of 210 min between orbital and rotational 
periods. Alas, for the present data, the fit is insensitive to a m , be- 
cause all exposures cluster around <p m = 0.25 or 0. 75. We fix dj n 
at 0.1 5, close to the pulse d fractions of blue l ight dHellier et al.l 
119871) and X-ray emission dRosenet al. 11199 ll) . 

Finally, we consider the emission properties of the irradiated 
atmosphere. The intensity emerging from an infinite isothermal 
layer with optical depth r along its normal is / oc 1 - exp(-r//j) 
with /j = cos 6. For large t, this form approaches the blackbody 
law and a realistic model should account for limb darkening or 
brightening. We find that all fits with free r yield r <; 10. Hence, 
we consider only the optically thick case that we describe by a 
linear intensity law / oc 1 - Mi.c a (l - /"), where u\ y ca is the limb 
darkening coefficient for the Call line that may be positive or 
negative. The contribution of surface element k to the line flux 
then is 

dF e . k = fl e ,l (iA k ^-^/j k (l - Ki,Ca(l -/*k))> (7) 

n 
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Fig. 10. Variation of the flux in the Call/18498 emissivity as a 
function of 77 = cos § for irradiation by a point source. The green 
line is for emission proportional to the energy influx, the black 
curve, for a simple analytical model that concentrates the emis- 
sion closer to the terminator, and the open circles for the numer- 
ical tomographic model. 



where a e ,ca,i is the proportionality factor that refers to the 
white dwarf as the irradiation source and the flux contribution 
is again set to zero for hidden elements (/^ = cosO^ < 0). 
Correspondingly, a e ,Ca,2 describes the 'hot spot' as irradiation 
source. The parameter cn m that shifts the limit of Call emission 
relative to the geometric terminator of each source is included by 
replacing j] k = cos # k in Eq. © by 77J. = (77k - ci im )/(0.9 - c iim ) 
that stays at unity for 77 = 0.9 (L; point) and vanishes at 77 = cum. 
We assume that the Nal emission shares all parameters with the 
Call emission except for its different normalization relative to 
Call expressed by a factor a e ,Na- The fit yields the parameters 
listed in Tab. [4] and, in addition, the values of the five propor- 
tionality factors a a , a ratio , a e ,Ca,i> fl e,Ca,2, and a e ,Na- 

6. Model fits 

Narrow emission lines have been detected in many C Vs and have 
been used to determine binary parameters assuming an origin 
from the irradiated face of the secondary. We confirm this ori- 
gin for EX Hya, but the interpretation of these lines is by no 
means straightforward. The centroid of the emission depends on 
the model parameters and minimum x 1 may occur for different 
values of the total mass M depending on the choice of f(rf) and 
cn m . The situation differs if absorption and emission lines are 
considered together because then the model covers the entire 
range of radial velocities that occur over the surface of the sec- 
ondary star. To obtain an internally consistent fit, we impose the 
side condition that the fit assumes minimum^ 2 for the Call lines 
and the Nal lines at the same M, by appropriately adjusting the 
parameters of the emission line model. All fit results presented 
in this paper comply with this condition. 

6.1. Ca/M8498 emissivity 

A representative set of emission models is shown in Fig.[l0] The 
simple assumption of emission proportional to the incident en- 
ergy flux, f{rj) oc 77 (green line), utterly fails yielding profiles at 
quadrature that lack intensity at the higher velocities and require 
relatively more emission at small 77. A one-parameter model of 




0.85 0.90 0.95 

M (M©) 



Fig. 11. Result of the multiparameter fit for the numerical model 
of the Call emissivities of Fig. [10] The Nal value (open cir- 
cles) is shown as a function of total mass M with the other pa- 
rameters free, the Call value Xc- d (+) is slaved to reach its min- 
imum at the same value of M (requiring c\\ m = -0.085). The 
green curve is a parabola fit and the dotted line indicates the 
99% confidence level. 



the form f(rf) oc 1 - exp(-ai]) fairs much better. The best fit 
requires a = 5.9 and cy lm = -0.060 (black solid line). The tomo- 
graphic model with eight fitted emissivities f m (Sect. [5]) yields 
a slightly improved fit. The displayed model (open circles) is 
for an intermediate level of regularization with a Lagrange mul- 
tiplier A — 100 and requires c\\ m = -0.085. Lower and higher 
values of A produce a more or less pronounced emission peak 
near 7/ 0.3 and the error bars on the data points reflect the 
range A = to 1000. Despite the slightly negative values of ci; m 
for both models, the emission is practically limited to the illu- 
minated part of the star if w e consider that the finite extent of 
the source at the white dwarf (Sie gel et al.| [l989) alone accounts 
for a fuzziness of the terminator of about 2° (Aci; m = ±0.035). 
Furthermore, the models lack the resolution to account for an 
abrupt drop of the emissivity at 77 = 0. Hence, the slight exten- 
sion of the emission beyond the geometric terminator is only 
marginally significant and we conclude that there is no com- 
pelling evidence for lateral energy transport in the atmosphere 
across the terminator. On the other hand, it is highly significant 
that the emission decreases much slower with decreasing rj than 
the incident energy flux. The observed increase in /(77V77 as the 
terminator is approached i s qu alitatively expected from the cal- 
culations of [Konig et al.J d2006l) for the hydrogen Lyman line 
emission from an irradiated (white dwarf) atmosphere. Realistic 
calculations are complex and need to consider the spectral en- 
ergy distribution of the incident radiation and its degradation as 
the terminator is approached in a spherically extended irradiated 
atmosphere treated in 2-D or 3-D. 

6.2. Line profiles and quality of the fit 

The high quality of the fit to the phase-resolved line profiles of 
the residual spectra (Sect. 13 - 31 > is demonstrated in Figs. [12] and 
[T3l Figures [T2h and d show excerpts of Fig. [5] that cover the 
NaL18183/8195 doublet and the CaIL18498 emission line, re- 
spectively, the model line profiles are shown in Figs. [TZb and 
e, and the residuals between data and model in Figs.[T2t and f. 
Figures [T3b,b display a quantitative version of the same result. 
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Fig. 12. a-f. Sections of Fig. [5] containing the NaI/18183/8195 doublet and the Call/18498 emission line. The Paschen background 
has been subtracted in the latter. The three panels each (a-c and d-f) show the data (left, a, d), the model (center, b, e), and the 
residual fluxes (right, c, f). The displayed intensities are relative to the orbital mean, the zero level is shaded gray, lower and higher 
intensities appear darker and lighter, respectively. FeI/18220.4 appears as very faint emission line in panels a and c, rows 7 to 10. 
The gray scales in the Nal and Call images differ. 
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Fig. 13. a,b. Comparison of observed phase-resolved line profiles (data points) and best-fit model spectra (green curves) for the 
NaI/18183/8195 doublet (left, a) and the Call/18498 emission line (right, b). The Paschen background has been subtracted in the 
latter. 



There is no background in the Call line profiles since they have 
been extracted from the profiles of Fig.[5]by removal of the un- 
derlying Paschen background. The enhanced emission at the ex- 
tremal velocities of the Nal doublet in Fig.[T2b are a result of the 
subtraction of the orbital mean spectrum from the data and the 
model as described in Sect. 13. 31 The sodium doublet is fitted be- 
tween 8160.5 and 8214. 5A and contributes 2880 data points in 
the 16 phase intervals, while the Call data set contains 250 non- 
zero data points. The best fit has^ a = 2608 and^ a = 274, or a 
total x 1 = 2882 for 3 1 16 d.o.f. With a reduced xl= 0.925, the fit 
is clearly good. The Nal part benefits from the inclusion of pres- 
sure broadening (Fig.[T3b), while the Call profiles (Fig.[T3b) are 



well matched without any additional broadening beyond radial 
velocity smearing. The humps that appear on the low-velocity 
slopes of the Call profiles at 098= 0.30 and 0.73 are responsible 
for the enhanced emissivity near rj = 0.3 in Fig. [10] Some vari- 
ability in the individual Call profiles can not be matched by the 
adjustment of parameters and may indicate temporal fluctuations 
of the line emission. 

Figure [14] shows the spectra observed from the illuminated 
and the dark side of the star shifted into their respective rest sys- 
tems. The emission line profile (top) displays dips flanking the 
emission peaks that represent residues of the incomplete fill-up 
of the underlying absorption components. The dips support our 
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Fig. 14. Observed (black) and fitted (green) line profiles of the 
Nal line for the illuminated (top) and dark side (bottom) of the 
secondary star. Same data as in Fig.Q~3]shifted into the respective 
rest system and co-added for two phase intervals (see Sect. 15. 4\ . 
See text for further explanation. 



choice of model for the composite Nal line profiles (Sect. 15. 2\ . 
The observed absorption line profile (bottom) is accompanied 
by model spectra with and without the orbital mean subtracted 
(solid and dashed green curves, respectively). They demonstrate 
that the subtracted model retains about 90% of the signal, pro- 
vided the spectral structure is restricted to a narrow wavelength 
range. The result justifies our use of the mean observed spectrum 
as a template for the telluric line correction. 

Finally, Fig. [7] compares the integrated observed and mod- 
eled Call emission line fluxes of Fig. [T3b as a functions of or- 
bital phase (upper panel, black histogram and green curve). The 
observed slight asymmetry of the light curve is modeled by as- 
suming the 'hot spot' as a second source of irradiation. Some 
other physical effect can not be excluded. 

6.3. System parameters 

Figure [TT] shows the variation of the %L vs. total mass M for 
the Nal lines with Xc a f° r the Call emission slaved to assume a 
minimum at the same M by appropriate choice of the parame- 
ters of the numerical emission model. There is a well defined^- 2 
minimum at a total system mass M = 0.898 M , with an error of 
0.026 M at the 99% confidence level (dotted line). This result is 
obtained with a radial velocity amplitude of the whi te dwarf of 
59kms _1 dBelle et al.ll2003l: iHoogerwerf et al.ll2004l) and a sys- 
tem inclination of 77.8°. The uncertainties in K\ and i of about 
3kms"' and 0.4° add +0.016 M G and +0.005 M to the error 
budget of M, respectively. Adding the errors quadratically, our 
measured system mass is M — 0.898 + 0.03 1 M . The remaining 
fit parameters are summarized in Tab. [4] 

The slightly mass-dependent inclination i = (77.8 ± 0.4)° 
is determined from the FWHM of the X-ray eclipse of 155 s 
(Muk ai et all 119981: IHoogerwerf etafl 12005b on the assump- 
tion that only the lower pole of the white dwarf is eclipsed 
and is responsible for the partial character of the X-ray eclipse 
(Beuermann & Osborne 1988; Beuermann et al. 2003). Some 
varian ts of the model dHellier et al.l 119871: iRosen et al~l 1 1988L 
[1991 may yield a slightly different value of i and we note that a 
±1° change in i causes a shift of +0.012 M in M. 



The velocity amplitude K 2 of the secondary star is a derived 
quantity and equals K 2 = (27rGM/P) 1/3 sin i - K\ = 432.4 ± 
4.8 km s _I . The mass ratio is also a derived quantity and given by 
q - K1/K2 = 0.1365 + 0.007 1 . The error in K 2 is almost entirely 
due to the noise in the line profiles, while the error in q reflects 
mostly the uncertainty in K\. The masses of primary and sec- 
ondary are M, = 0.790 ±0.026 M andM 2 = 0.108 ±0.008 M , 
respectively. The mean Roche lobe radius of the secondary star is 
R 2 = 0.1516 ± 0.0034 R . Our direct fit to the line profiles elim- 
inates all problems associated with the determination of radial 
velocities from the complex line profiles as far as possible. Our 
result for K 2 significantly exceeds t he value of 360 ± 35 kms re- 
ported bv lVande Putte et aT71 d2003l) . on which most recent pub- 
lished Mi values were based. The system velocity obtained from 
the fit is y — -58.2 ± 1 .0km s _1 in agreement with the result ob- 
tained in Sect. |4.2| from the radial velocity curves. 

6.4. Roche tomography 

Roche tomography subjects the emission of all surface elements 
to some type of regu larization and works bes t if the binary pa- 
rameters are known dWatson & Dhillonll2004l) . Since the deter- 
mination of these parameters is our main goal, our tomographic 
approach involves the following simplifications: (i) there is no 
freedom in the contributions of the individual surface elements 
to the Nal absorption line profiles (Sect. 15.31 ); and (ii) the con- 
tributions to the emission line profiles depend only on the angle 
of incidence 1? of the irradiation, with shadowing by the accre- 
tion disk superimposed (Sect. |5.4T i. However, as discussed above, 
there is freedom in the dependence of the emission line flux on 
the incident energy flux, in the extent to which the emission ex- 
tends into the 'dark' side of the star, and in the angular distri- 
bution of the emission line intensity that emerges from a given 
surface element (Sect. [5T4b . The fact that the model fits the data 
with a reduced ;^=0. 92 indicates that any added freedom in the 
tomography can not be expected to improve the fit. Such an ap- 
proach would need data of much better statistical significance. 

Figure [15] shows the images of the secondary star derived 
from our restricted tomography as seen at orbital phase 0.75 and 
an inclination of 77.8°. The left panel depicts the Nal image and 
the right panel the Call image. The pictures are based on the nu- 
merical model of the emission line fluxes of Fig. [10] The princi- 
pal features are the clear distinction between dark and irradiated 
hemispheres of the secondary star and the dark lane representing 
the shadow cast by the accretion disk. The geometric terminator 
for illumination by a point source at the position of the white 
dwarf is indicated by the dashed line. Emission is seen to cease 
rather abruptly within a few degrees of the terminator, except for 
the faint emission produced by the ' hot spot' that extends beyond 
the terminator and fills in the disk shadow. In the Nal image, the 
effect of limb (and gravity) darkening is visible as a reduced ab- 
sorption line flux as the limb is approached. The shadow of the 
accretion disk takes away emission very close to the Li point, 
but otherwise has little effect on the relative distribution of the 
emission. 

7. Discussion 

The main result of our high-resolution spectrophotometric study 
is the discovery of narrow absorption and emission lines from the 
photosphere and chromosphere of the illuminated secondary star 
in EX Hya. We present a novel method for the kinematic analysis 
that involves a direct fit of a model of the emitting Roche-lobe 
filling secondary star to the observed absorption and emission 
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Fig. 15. a,b. Tomographic images of the secondary star of EX Hya seen at <pgs = 0.75 and i = 77.8° for the numerical model 
of the emissivities of Fig. [10] (a) Left: NaI/18183/8195 absorption/emission, (b) right: Call/18498 emission. The x' - z' plane is 
perpendicular to the line of sight and coordinates are relative to the binary separation a with the origin in the white dwarf. The 
geometric terminator for illumination by a point source at the white dwarf is indicated by the dashed line. The intensity scales differ 
for both images. In the left hand panel, black is absorption, white emission, and the red rim represents zero intensity, while the 
dark side of the star in the right panel denotes zero line flux that increases over red to white. Faint emission in the disk shadow and 
beyond the terminator is due to irradiation by the 'hot spot'. 



line profiles. A key quantity is the linear size of the Roche lobe 
between its back and the Li -point (Sect. l4~3l l. which we deter- 
mine from the combined fit of the synthesized line profiles to 
the absorption and the emission line data. Combined with the 
previousl y measured radial velocity amplitude K \ of the white 
dwarf dBelle et al.ll2003b iHoogerwerf et alj|2004l) . this line syn- 
thesis approach allows us to determine accurate masses of the 
binary components. 



7. 1 . The secondary star 

Our fit yields a secondary mass M 2 = 0.108 ± 0.008 M Q and 
a mean Roche radius R 2 = O.1516R =1.O5 x 10 10 cm, which 
place the star much closer t o the main sequence tha n the former 
estimate of 0.078 M G (e.g., iBeuermann et alj|2003l) . The abso- 
lute #-band magnitude of EX Hya B is M K = 8.83 + 0.12 in 
perfect agreement with the spectral classification dM5.5 + 0.5. 
The derived parameters of EX Hya B are summarized in Tab. |5j 
with the data of the nearby dM5. 5 star G1551 that has an in- 
terferometrically measured radius (Segransan et al. 2003) added 



Table 5. Parameters of the secondary star in EX Hya compared 
with those of G1551 (Proxima Cen). 



Parameter 


EX Hya B 


G1551 




Spectral type 


dM5.5 


±0.5 


dM5.5 




Temperature T tfl (K) 






3042=t 


117 


Distance d (pc) 


64.5 ± 


1.2 


1.295 : 


t 0.007 


K (mag) 


12.89 = 


= 0.11 


4.36 




M K (mag) 


8.83 ± 


0.12 


8.80 ± 


0.02 


Mass M (M Q ) 


0.108 = 


= 0.008 


0.123 : 


t 0.006 


Radius R (Rq) 


0.1516 


± 0.0034 


0.145 : 


t 0.011 


Velocity amplitude K 2 (kms~') 


432.4 = 


= 4.8 






Systemic velocity y (km s~' ) 


-58.2 


t 1.0 







for comparison. The radius of an 0.1 08M H model star o f so- 
lar composition is 0.925 x 10 10 cm dBaraffe et al. 1 1 1998b and 
rise s to 0.98 X 10 10 cm if rotational deformation is accounted 
for dRenvoize et al. Il2~002l) . These models apply to stars without 
spots, while spotted stars tend to have somewhat larger radii. 
Just how much this effect influences the radii of the secondary 
stars in CVs is not known (see Beuermann et al. 2006 and ref- 
erenc es therein). The calibr ation of the radii of field stars vs. 
M K (IBeuermann etal.1 1 19991 Eq. 7) yields R = O.1514R for 
Mk = 8.83, confirming that EX Hya B is essentially a main 
sequence star, but differs from a field star of the same absolute 
magnitude by being deformed and irradiated. 

7.2. The white dwarf 

Previous estimates of the masses of the two stellar compo- 
nents in EX Hya relie d on the uncertain veloci ty amplitude of 
the secondary star of IVande Putte et al. I (l2003h . K 2 = 360 + 
35kms _I . Our result, K 2 = 432.4 + 4.8 km s -1 , is 2.1cr from 
Vande Putte's result and yields a much more accurate kinematic 
solution. The primary mass Mi = 0.790 + 0.026 M super- 
sedes the value of M\ ^ 0.46 M quoted in most recent stud- 
ies on EX Hya (e.g. Beuermann et al] 2003 : Vande Putte et al. I 



l2003r.lHoogerwerf et alJ2004l2005l:lMhlahlo et al. l2006h l'L The 
new mass is close to the mean for short-period CVs, (M) ^ 
0.8 M Q , and confirms the early result M\ = 0.78 ± 0.17M Q 
of iHellier et alJ (Il987l) . Given the larger primary mass, one 
ca n re-estimate the mass transfer rate following the analysis 
of IBeuermann et all d2003l see also Ritter (1985)). From their 
Fig. 1, M\ is seen to drop with increasing Mi and reach about 
the value expected from gravitational radiation as the sole mo- 
mentum transfer process for Mi = 0.79 M . 

The radius of a white dwarf of 0.790 + 0.026 M with an in- 
trinsic effective temperature of about 15000 K is Ri = (7.35 ± 



From Segransan et al. (2003) 



1 The primary mass of 0.91M o derived by iBelle et al.l {2003) from 
Ki = 59.6 kms~', K 2 = 360 kms~', i = 78°, and an assumed M 2 = 
0. 152 M violates Kepler's law. 
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0.23) x 10 8 cm ba sed on models with a thick hydrogen enve- 
lope (Wood 1995). The mean temperature of the white dwarf 
in EX Hya determined from the HST FOS spectrum is 25000 K 
(Eisenbart et al. 2002), but that temperature is dominated by the 
heated polar caps responsible for the pronounced spin modula- 
tion seen in X-rays and in the UV and the temperature of the 
underlying white dwarf is likely to be lower, roughly as noted 
above. With the quoted radius, the gravitational redshift at the 
surface of the white dwarf is y grav = 47.6 ± 3.1kms _1 . The 
predicted apparent systemic velocity of the white dwarf then is 
7 + fgrav = ~10.6 + 4.1 kms _I . The velocity derived from the 
X-ray emission lines that origin ate close to the white dw arf sur- 
face is y — -2.8 ± 2.3 km s -1 dHoogerwerf et al.ll2004l) . where 
the error is the statistical one and the systematic error is larger 
(C. Mauche, private communication). While the interpretation of 
these numbers in terms of a gravitational redshift measurement 
may be premature, it is clear that an independent measurement 
of the white dwarf mass becomes feasible with a more secure 
value of the apparent systemic velocity of the white dwarf. 



Watson, C. A., & Dhillon, V. S. 2004, AN, 325, 189 
Wood, M. 1995, LNP 443, 41 
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